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We show that single-component fermionic polar molecules confined to a 2D geometry and dressed by a mi- 
crowave field, may acquire an attractive 1 /r 3 dipole-dipole interaction leading to superfluid p-wave pairing at 
sufficiently low temperatures even in the BCS regime. The emerging state is the topological p x + ip y phase 
promising for topologically protected quantum information processing. The main decay channel is via col- 
lisional transitions to dressed states with lower energies and is rather slow, setting a lifetime of the order of 
seconds at 2D densities ~ 10 8 cm -2 . 



Remarkable progress in the studies of ultracold atomic 
Fermi gases (HQ] has opened up prospects for creating novel 
phases of fermionic atoms. Of particular interest is the topo- 
logical superfluid p x + ip y phase for identical fermions in two 
dimensions (2D) [3], discussed in the contexts of superfluid 
3 He and the fractional Quantum Hall effect 0, 0. The in- 
tense interest arises from the exotic topological properties of 
the phase at positive chemical potential /i > (i.e. in the BCS 
regime). In the presence of vortices, the groundstate becomes 
highly degenerate, spanned by zero-energy Majorana modes 
on the vortex cores O.^]. The highly non-local character of 
these states is expected to suppress decoherence processes, 
and allow this degenerate subspace to be used for topologi- 
cally protected quantum information processing 0]. 

The p x + ip y topological phase has been predicted to be 
the groundstate of ultracold fermionic atoms interacting via 
a p-wave Feshbach resonance [3]. However, the realization 
of the p x + ip y phase in this way encounters serious difficul- 
ties. Away from a Feshbach resonance the superfluid tran- 
sition temperature is vanishingly low. While it may be in- 
creased on approach to the resonance, in this case the system 
becomes collisionally unstable. Fermions form long-lived di- 
atomic quasibound states and their collisions with the atoms 
cause relaxation into deep molecular states, leading to a rapid 
decay of the gas SH. 

In this paper we show that a stable topological p x + ip y 
phase can be created with fermionic polar molecules with 
large dipole moment. Ultracold clouds of polar molecules 
in the ground ro-vibrational state have been obtained in re- 
cent successful experiments Jiol [Till . Fermionic 40 K 87 Rb 
molecules [10] have a permanent dipole moment d ~ 0.6 
D, and the dipole moment of 6 Li 133 Cs fermionic molecules 
should be close to 6 D, the same as for the created bosonic 
molecules 7 Li 133 Cs 111 ill . Being electrically polarized such 
molecules interact via long-range anisotropic dipole-dipole 
forces, which has crucial consequences for the nature of quan- 
tum degenerate regimes. In particular, this provides the pos- 
sibility of superfluid pairing at sufficiently low temperatures 
in a single-component Fermi gas. In 3D the ground state 
of a gas of fermions with dipole moments aligned in the z- 
direction has a pairing function which vanishes for p z = 
1 12I1 . The 2D Fermi gas of canted dipoles has a ground state 



with a pairing symmetry of a similar form 111 311 . In both 
cases, the presence of nodes in the order parameter make these 
phases distinct from the p x + ip y topological phase. Our route 
to a stable p x + ip y phase is somewhat simpler than other 
approaches [lj]. 



Our idea is to use polar molecules confined to a 2D ge- 
ometry and dressed by a microwave (MW) field which is 
nearly resonant with the transition between the lowest and the 
first excited rotational molecular levels. As we describe be- 
low, the dressed polar molecules acquire an attractive l/V 3 
dipole-dipole interaction, which leads to superfluid pairing of 
Px + ipy symmetry. Staying in the BCS limit, the superfluid 
transition temperature can be made sufficiently large, and de- 
cay processes sufficiently slow, to allow realization of this 
phase in experiment. The effects of MW dressing of polar 
molecules have been considered as a way to tune the inter- 
molecular potential 11511 . and to form a repulsive shield for 
suppressing inelastic losses lflrjll . The possibility of attractive 
interactions, which we consider here, has very important con- 
sequences for the nature of the phases that can arise, and for 
the stability. 

We consider a gas of fermionic polar molecules which are 
tightly confined in one (z) direction and assume that the con- 
finement length l z still greatly exceeds the size of a molecule. 
Then the translational motion of the molecules is 2D, but rota- 
tional eigenstates \J, Mj) are those of the 3D molecule. The 
operator of the dipole moment d can have non-zero matrix el- 
ements only between states with different rotational quantum 
numbers J. The transition dipole moment for J = — > J = 
1 is dt = | (0, 0|d|l, Mj) | =d/y/3, with Mj = 0, ±1, and d the 
permanent dipole moment of the molecule. 

We then apply a circularly polarized MW field which prop- 
agates in the z direction and has a frequency uj close to the 
frequency oj of the transition between the states |0, 0) and 
|1,1). If the Rabi frequency fi^ = d t E/T% and the detun- 
ing 5 = uj — ujo satisfy the inequality |5|,f2fl <C u>o< then 
the rotating wave approximation is valid and the MW electric 
field E[t) couples only the states |0, 0) and |1, 1). The result- 
ing states may be represented in the dressed-molecule picture, 



2 



with wavefunctions lfl7l 

|+) = a\0,0;N) + be- iu,t \l,l;N-l) 



0.125 



= b\0,0;N) 



ae- luJt \lA:N -1) 



(1) 

(2) 



where N labels the number of photons in the field, and 
a = -A/^A 2 + fi£, b = Vl R /^/VL z R + A 2 , A = (S + 
yjb 2 + ifl 2 R ) /2. We will consider S > ^b. , and choose S > 
such that the energy of the state |+) lies above the energies of 
|— ) and |1, —1). If the MW field is ramped on adiabatically, 
then the ground state |0, 0) evolves into the state |+), and all 
molecules can be prepared in this state. As described below, 
relaxation to the lower energy states, |— ) and |1, — 1), can be 
very slow. 

We derive the interaction potential between two molecules 
within the Born-Oppenheimer approximation, in which the 
molecules are assumed to be at fixed locations with a sep- 
aration r = r(cosc/), sm<j)). At large separations, the 
molecules are both in the state |+). Each molecule has an 
effective electric dipole moment d c g = — \f2abd t , which 
rotates in the plane of translational motion: (+|d|+) = 
d e s(cosu>t, smut, 0). The interaction potential at large dis- 
tances is then 

nO = ^^^^ = %[l-3cos 2 M-0)]. (3) 



Thus, the time-averaged interaction is attractive: 

Fo(r-oo) = -<4/2r 3 



(4) 



and is characterized by the lengthscale r* = Md 2 s /2T% 2 , 
where M is the mass of a molecule. The quantity r* is de- 
fined analogously to the van der Waals length for atoms and is 
a measure of the radius of the centrifugal barrier experienced 
by the (fermionic) molecules. 

At smaller separations, the dipolar interactions between the 
molecules cause them to depart from the state |+). This 
occurs when the characteristic interaction energy d 2 / r 3 be- 
comes larger than the detuning %\5\, setting a new length- 
scale rg = [d 2 /(HIS])] 1 / 3 . We have found the resulting Born- 
Oppenheimer surfaces using a full coupled channel calcula- 
tion containing all the levels | J, Mj; N). Similar calculations 
are described in Refs. Il5l [lal . We choose a positive detun- 
ing 5 > and assume that the lengthscale of the potential, 
rs, is larger than the confinement length l z , so that the inter- 
action is 2D. The potential energy curves of even parity are 
illustrated in FigQ] showing a potential Vo(r) that has a re- 
pulsive core for r < rs and is attractive at r > rs, with a 
long-range 1/r 3 tail. As discussed below, the repulsive core 
prevents low-energy particles from approaching each other at 
distances r < rs and suppresses inelastic collisions, including 
"ultracold chemical reactions" recently observed at JILA for 
KRb molecules. 

We now analyze the low-temperature phase of a 2D gas of 
identical fermions interacting via the potential Vo(r), assum- 
ing that rs <C r* . Due to the presence of an attractive 1/r 3 
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FIG. 1: Potential energy curve Vo(r) for two |+) state molecules, 
computed for fin. — 0.255 (see text). Anticrossings with other field- 
dressed levels of even parity occur at distances r ~ rg , as shown in 
the inset. 



tail given by Eq. |@), one expects that the Fermi gas is unstable 
to the formation of a superfluid state. In the ultracold dilute 
limit, where the momenta of colliding fermionic particles sat- 
isfy the inequality kr* -C 1, this tail provides a contribution 
oc (—kr*) to the scattering amplitude. This is the so-called 
anomalous contribution coming from distances of the order 
of the de Broglie wavelength of the particles and obtained in 
the Born approximation I19I1 . It greatly exceeds the leading 
short-range (r < r*) contribution which is related to the p- 
wave scattering and is oc k 2 away from p-wave resonances. 
Thus, omitting second order corrections, a detailed behaviour 
of the potential Vo(r) at distances r < r* drops out and the 
only important lengthscale is r* . 

In our analysis of the superfluid phase we confine ourselves 
to the BCS weak coupling regime, where kpr* <§; 1 with 
kp = \fAim being the Fermi momentum and n the gas den- 
sity. (We consider a uniform 2D gas; effects of a trap can be 
included within the local density approximation.) The regular- 
ized gap equation is obtained expressing the interaction poten- 
tial through the zero energy vertex function r(fc, q) governed 
byH 



T(k,q) = V (k-q)-J 0, 



d 2 q' T(k,q')V (q-q') 



E a 



(5) 



with E q = h 2 q 2 /2M, and Vo(q) being the Fourier transform 
of Vo(r). The gap equation then reads J2, 2(1 21]: 



where 



(2n) 2l[k ' q) 2 



tanh(e q /2T) 



E a 



(6) 



J (E q — /i) 2 + | A q | 2 is the energy of single- 



particle excitations, and \i > is the chemical potential which 
is equal to the Fermi energy Ep — h 2 k F /2M. To first order 
we replace T(k, q) in Eq. © with V (k — q). At T — we 
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put tanh(e q /2r) = 1 and perform an analytical analysis as- 
suming that in the weak coupling limit the main contribution 
to the integral in Eq. (0 comes from momenta q close to kp. 
It shows that the dominant pairing instability is in the channel 
with orbital angular momentum 1 = 1. The most stable low 
temperature phase has p x ± ip y symmetry, following from the 
fact that this phase fully gaps the Fermi surface, in contrast to 
competing phases lEill . A full numerical solution of the reg- 
ularized gap equation confirms this analysis. It further shows 
that | Afe| rises linearly for fc < fc F , and approaches a constant 
— E F exp(-37r/4fcf r*) for fc > kp. 

In the 2D geometry that we consider, the critical tempera- 
ture T c of a Fermi gas is set by the Kosterlitz-Thouless tran- 
sition. However, in the weak coupling limit the Kosterlitz- 
Thouless temperature is very close to T c obtained in the BCS 
approach B23I1 . For T — > T c we omit |A q | in the expression 
for e q in Eq. © and obtain: 



T c w £ F exp(-37r/4fc F r*), 



(7) 



where the numerical prefactor is of order of unity [24]. Thus, 
to obtain an achievable value of T c one requires k F r* to be 
not much smaller than unity. The BCS approach assumes that 
the exponential factor in Eq. © is small and T c <C Ep. A 
limitation on the strength of the attractive interaction is set 
by the condition of stability to phase separation (collapse to 
a high-density gas). A full calculation of this limit requires 
a strong-coupling theory. However, estimates (provided by 
Hartree-Fock theory) suggest that the compressibility is pos- 
itive for k F r* < 3.7. Thus, the regime of moderately strong 
interactions k F r* ~ 1 is accessible. 

The p x + ip y phase spontaneously breaks time-reversal in- 
variance (the phase p x — ip y is its degenerate time-reversed 
partner)[25]. It can be viewed as a state in which the Cooper 
pairs have an orbital angular momentum of Ti with respect to 
the z-axis. The p x + ip y phase can exist in one of two topo- 
logically distinct phases, depending on the sign of the chem- 
ical potential[4]. The phase at [i < may be continuously 
deformed to the vacuum state; the phase at u > is topologi- 
cally distinct from the vacuum and has several very interesting 
properties. Most notably, the vortices of this phase carry local- 
ized zero-energy states, described by a Majorana fermion on 
each vortex core. These lead to non-abelian exchange statis- 
tics ASlal and possible applications for topologically protected 
quantum information processing[7]. The superfluid of dipo- 
lar interacting spinless fermions that we describe above has 
/i > and is in the relevant topological phase. 

The typical interatomic potential between atoms or 
molecules (without MW dressing field) has a short range 
i?o ~ 1 — 10 nm. Thus, in 2D the scattering phase shift 
is ~ (k F R ) 2 so that T c ~ Ep exp[— 1/ (kpR a ) 2 ] is van- 
ishingly small. If the interaction strength is tuned close to a 
Feshbach resonance [3], such that the transition temperature 
strongly increases, then the particles have a significant proba- 
bility to be inside the centrifugal barrier at separations of the 
order of Rq. Under these conditions the system is very sus- 
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FIG. 2: Inelastic rate constant a as a function of r$ /\s for fi^ 
0.255 and kr* = 1 (see text). 



ceptible to rapid losses arising from collisional relaxation into 
deep bound states|H,|3l. 

In contrast, for MW dressed molecules interacting via the 
potential Vo(r), the contribution kpr* to the attractive cou- 
pling strength provides a significant transition temperature T c 
even far from the resonance associated with the presence of 
a two-molecule bound state. Given that the molecules have a 
small probability to be at separations ~ rg, in this BCS regime 
one anticipates that the superfluid phase is not susceptible to 
relaxation losses. 

The dominant loss mechanism is from binary inelastic 
collisions between |+) molecules, in which one or both are 
transferred to the state |— ) or |1, — 1), which (since 5 > 0) lie 
lower in energy than |+). For fi F < S, the released kinetic 
energy is ~ %8 and can cause both molecules to escape from 
the sample. The kinetic energy release requires a momentum 
transfer of ~ Ti/\$ with \$ = y/h/M5. For \g/r$ <C 1 the 
particles cannot approach each other sufficiently closely to 
allow the required momentum exchange, and one anticipates 
a reduction in the loss rate. The same condition can be derived 
semiclassically as the condition of adiabatic motion in the po- 
tential. To go beyond this limit, and determine the loss rate for 
\s ~ rs, we have solved the full two-body scattering problem, 
involving states of even parity which at infinite separation are: 
(|+>,|+», (|+>,|-», (+>,J1,-1», (|->,|-», (|->,|1,-1» 
[the state (|1, — 1), |1, — 1)) is decoupled]. We calculate (nu- 
merically) the probabilities Pi that two |+) -state molecules 
with relative angular momentum / are scattered into any 
outgoing channel in which at least one of them is in the state 
|— ) or |1, —1). This corresponds to non-adiabatic transitions 
from the potential Vq (r) to the other potentials shown in the 
inset to FigfT1ll26ll 

Taking into account that two molecules are lost in each in- 
elastic collision, and writing the molecule loss rate as n = 
—an 2 , for the 2D inelastic rate constant we obtain: a = 
4h/M J^i Pi- F° r incident energy chosen such that kR* <C 1, 
we have a oc (kr*) 2 . The results for this quantity versus 
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rs/Xs are shown in FigfS] for Q/S = 0.25, and kr* = 1. The 
general trend is a reduction of inelastic losses with increasing 
rs/Xs, consistent with the semiclassical expectations. How- 
ever, in addition there is a dramatic modulation of the in- 
elastic scattering rate, arising from an interference of incom- 
ing and outgoing waves in the scattering potential. By tun- 
ing to rs/Xs ~ 10.5, the rate constant can be suppressed to 
a~4x IQ-^Ti/M. 

Thus, at a density n = (10 8 — 10 9 )cm~ 2 of, for example, 
7 Li 40 K molecules the lifetime of the gas is r i=s (an)^ 1 ~ 2 — 
0.2 s. The permanent dipole moment of 7 Li 40 K in the ground 
state was found to be 3.5 D 02711 . and for rs/Xs ~ 10.5 and 
fin = 0.255 the lengthscales are r$ ~ 30 nm, r* ~ 200 nm. 
For kpr* close to unity we then get n~2x 10 8 cm -2 and 
Ep ~ 120 nK, so that the transition temperature is T c ~ 10 
nK and the lifetime is ~ 1 s. (For 40 K 87 Rb, the possible r* is 
rather small and the high densities required for a sizeable T c 
lead to rapid losses. The addition of a shallow optical lattice 
will increase the effective mass M, allowing T c ~ 10 nK at 
n ~ 10 8 cm -2 .) 

We should avoid the presence of bound states of two 
molecules in the potential Vq (r), otherwise three-body recom- 
bination will lead to a rapid decay of the gas on approach to 
the superfiuid transition. A dimensional estimate for the three- 
body decay rate gives r" 1 ~ (for* 2 /M)(kpr*) 4 n 2 , which 
can be large for kpr* ~ 1 and reasonable densities. However, 
for S}r ~ 0.25(5 (as used above) the potential Vo(r) does not 
support bound states for rg/ X$ < 14. Thus, for the considered 
value rs/Xs — 10.5 the three-body recombination is absent. 

The formation of the p x + ip y superfiuid phase should be 
apparent in numerous observables. These include quantities 
used to detect s-wave pairing in two component Fermi gases, 
such as the density distribution, collective modes, and RF ab- 
sorption spectra Illy]. The most striking new features of the 
Px + ipy superfiuid arise in the presence of quantized vortices, 
which may be generated by rotation of the gas as in usual 
superfiuids. RF absorption will then show evidence for Ma- 
jorana modes on the vortex cores 12811 . Ultimately one would 
hope to probe non-abelian exchange statistics of these vortices 
& 
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